Spatial Inhomogeneities in Disordered d-Wave Superconductors: Effect on Density of 

States and Superfluid Stiffness 
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We study a short coherence length d-wave superconductor with finite density of unitary scatterers using 
the Bogoliubov-deGennes technique. We find that the low-energy density of states is reduced, the superfluid 
stiffness is significantly larger and off-diagonal long range order is more robust than the self-consistent T- 
matrix prediction. These results are a consequence of the inhomogeneous pairing amplitude in the ground 
state and of the low-lying excitations formed by hybridized impurity resonances. These features, with their 
nontrivial spatial structure, cannot be adequately described within the conventional T-matrix approach. 
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Recently there has been a great deal of theoretical 
and experimental interest in the problem of disordered 
d-wave superconductors (SC) given the d-wave symme- 
try of the high Tc cuprates. A single unitary scatterer 
was predicted by Balatsky and coworkers [|] to lead to a 
low-energy resonance with a characteristic four-fold sym- 
metric wavefunction about the impurity site. This was 
recently observed in an STM study U of a Zn-doped 
cuprate. The T-matrix approximation used in Ref. |^] 
is very accurate for the one-impurity problem, and the 
order parameter suppression near the impurity, which it 
neglects, does not lead to any qualitative changes. On 
the other hand, the problem of a finite density of uni- 
tary scatterers is more subtle There is a large body 
of theoretical work using the self-consistent T-matrix ap- 
proximation leading to very interesting predictions [^-0 . 
However, the impurity averaging procedure used in this 
method, and in other approaches 1^ that go beyond it, 
basically washes out the inhomogeneous structures that 
the disorder potential gives rise to. 

We use the Bogoliubov-deGennes (BdG) approach to 
study the spatial inhomogeneity induced by unitary scat- 
terers in a short coherence length (^o) superconductor, 
and to understand how these affect the low energy prop- 
erties of the system such as the one-particle density of 
states (DOS) N{uj) and the superfluid stiffness D^- Our 
main results can be summarized as follows: 

(1) The low energy DOS is considerably reduced relative 
to the T-matrix result. The low lying excitations are 
found to be generated by the interference of individual 
impurity resonances. Such excitations, with their non- 
trivial spatial structure, cannot be adequately described 
within the T-matrix formulation. 

(2) The superfluid stiffness Dg is found to be significantly 
larger than that obtained within the T-matrix analysis. 
The larger Dg directly correlates with the lower DOS for 
"normal" excitations. 

(3) We find that off-diagonal long range order (ODLRO) 
and finite superfluid stiffness survive to impurity concen- 
trations much higher than the critical concentration of 



the T-matrix approximation. This relative insensitivity 
of short coherence length d-wave superconductors to im- 
purities is shown to be closely tied to the inhomogeneity 
of the pairing amplitude on the scale of in response to 
a random potential, fn contrast, the T-matrix approach 
assumes a uniform amplitude which then gets globally 
suppressed to zero at a critical disorder. 

Several authors have previously used the BdG ap- 
proach for dirty d-wave systems. Tc reduction, super- 
fluid density and localization of excitations was stud- 
ied in Ref. |^ and more recently the density of states 
has been studied in refs. [ p^|jTl| ]. Our calculations dif- 
fer from these in several aspects, some more important 
(e.g., working at fixed density rather than fixed chem- 
ical potential and inclusion of inhomogeneous Hartree- 
Fock shifts) than others (such as choice of Hamiltonian, 
parameters, and particle- hole asymmetry). While our 
results are broadly consistent with those obtained previ- 
ously, what is new here is our emphasis on understanding 
the BdG results for N{uj), Ds and ODLRO in terms of 
two different effects: (A) the inhomogeneity in the lo- 
cal pairing amplitude which characterizes the disordered 
ground state, and (B) the spatial structures characteriz- 
ing the low-lying excitations in the disordered system. 
This provides a deeper insight into our BdG results, 
and also highlights the shortcomings of the T-matrix ap- 
proach. 

We model the 2D disordered d-wave SC by the Hamil- 
tonian Ti = IC + Hint + Tidis- The kinetic energy JC = 
—tJ2<cij> ai'^ia'^ja + h.c.) describcs electrons, with spin 
a at site i created by c]^, hopping between nearest- 
neighbors < J/j > on a square lattice. The inter- 
action term [|2| Hint = JJ2<ij> (^i ■ Sj - n^rij/A) + 
^X^i^iT^U is chosen to lead to a d-wave SC ground 
state in the disorder-free system. The spin operator 
— ^ia'^ap^iP^ whcrc the (7° are Pauli matrices, and 
the density riia — cla'^ia with rii — n.i| 4- Jiij. Finally, 
Hdis — X^i (^(*) ~ m) where ^ is the chemical poten- 
tial and the disorder potential V{i) is an independent 



1 



random variable at each site which is either +Vb, with a 
probabihty nimp (impurity concentration), or zero. We 
beHeve that such a simple model is adequate to describe 
the strongly-correlated cuprates at low temperatures be- 
cause their SC state has sharp quasiparticle excitations. 
The BdG equations are given by: 

where |u„(.?) = -J2sit+Wj)u„{j + 5) + iVij)-tl,)u„ij) 
and Au„(j) = J2s ^0 ~^ + similarly for 

v„{j). The pairing amplitude on a bond {j',5), where 
S = ±x, ±y, is defined by A{j;S) = —J{cj+siCj-\ + 
CjiCj+si) /2. The inhomogeneous Hartee-Fock shifts are 
given by flj = n - U{nj)/2 -\- JY^giuj+s) and Wj = 

We numerically solve for the BdG eigenvalues 
En > and eigenvectors on a lattice of 

N sites with periodic boundary conditions. We 
then calculate the pairing amplitude A(j;5) — 
JEn K(J + + ^nUKU + S)] /2 at T = 0, the 

density (hj) — 2^jjw„(j)p, and Fock shift Wj — 
JJ2n ^"(j + ^)^nij)- These are fed back into the BdG 
equation, and the process iterated until self consistency 
|l4| is achieved for each of the (local) variables defined on 
the sites and bonds of the lattice. The chemical poten- 
tial n is chosen to obtain a given average density (n) — 
Yi{ni) /N , The d-wave pairing amplitude is given by 
A(j) = [A(j; +x) - A{j; +y) + A(j; -x) - A(j; -y)] /4. 

We have studied the model for a range of parameters 
and lattice sizes. Here we focus on J = [/ = 1.15, in units 
of t = 1, with (n) — 0.875 (similar to the parameters 
used in refs. P,p^) on systems of size up to 26 x 26. For 
these parameters, and rii^p — 0, the DOS A^o — 0.21 and 
Ao ~ 0.077 corresponding to a maximum gap of 0.31. 
For the impurity potential we choose Vq = 100, close to 
the unitary limit. The results are averaged over 15-40 
different realizations of the random potential. 

Let us first study the density of states (DOS) 
Niio) = [\u,M'S{u;-En) + \vn{i)mu; + En)] 

(where we broaden the delta functions with a width com- 
parable to average level spacing). In Fig. 1 we plot N{lu) 
for several impurity concentrations on a small energy 
scale; for comparison, the maximum energy gap in the 
disorder-free system is 0.31 and the T-matrix self energy 
scale 7= ^ni,„pA/2A^o < 0.25 for the parameters 
chosen; (A is the T-matrix gap) . In the T-matrix theory 
N(lu) is a constant for w < 7, while we find a sharp dip 
in the DOS close to the chemical potential, consistent 
with rcf. In fact, we found A^(0) — for each impu- 
rity configuration at every concentration that we studied. 
The scale of the sharp dip at finite nimp was found to be 
the same as the energy of an isolated impurity resonance. 

It is very clear that the low energy DOS in the BdG 
calculations is considerably smaller than that in the T- 
matrix approximation (even though we do not have the 
spectral resolution to quantify the asymptotic form of 



N{uj) as w ^ 0). To highlight this, we compare in Fig. 2 
the finite 7V(0) of the T-matrix analysis ||l|| with the BdG 
N{0), which is the average of N{uj) over the (arbitrarily 
chosen) range < 0.05 <^ 7. 
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FIG. 1. Density of states (DOS) on a iV = 24 x 24 sys- 
tem, with J = U = l.lbt and (n) = 0.875, averaged over 40 
disorder realizations at each riimp. Note the sharp drop in the 
DOS near a; = on a scale much smaller than the energy gap 
of O.Slt in the pure system. 
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FIG. 2. BdG density of states (DOS) iV(0), defined as the 
average of N{u;) over the range < 0.05, is much smaller 
than the corresponding T-matrix result. The parameters are 
the same as in Fig. 1 and the normalizing factor is the pure 
system DOS iVo = 0.21. 

To gain further insight into this difference between the 
T-matrix and BdG results, we study the wavefunctions 
of the low-lying excitations for individual disorder real- 
izations. The probability density |w„(«)p -I- |'!;„(i)p corre- 
sponding to the lowest energy states at various impurity 
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concentrations are plotted in right hand panels of Fig. 3. 
Single impurity 
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FIG. 3. Left column: Evolution of the local pairing ampli- 
tude A(i) with impurity concentration. Dark regions in the 
grey-scale plot indicate suppressed pairing amplitude, and are 
correlated with the impurity locations. Parameters used are 
J = U = 1.15t and (n) = 0.875 on an = 24 x 24 system. 
Right column: The corresponding probability density 
|'itn(4)P + |i'n(i)|^ for the lowest excited state (n — 1) wave- 
function. Higher probability is indicated by a darker shade. 
Each impurity location is marked by a dot. 

The resonance for a single unitary impurity shows char- 
acteristic powerlaw tails along diagonal directions ||l|,|l6| . 
From Fig. 3, and other low lying excitations not shown 
here, we see that for finite nimp these wave functions 
are generated by the hybridization of individual impurity 
resonances. The effects of constructive and destructive 
interference between the "diagonal tails" of individual 



resonances are apparent. The importance of such states 
was suggested in ref. however, their analysis assumed 
that the resonance energies are randomly distributed over 
a scale W ^ Ao, which is not the case in the physical 
situation obtained here. 

We emphasize that excitations with such non-trivial 
spatial structures cannot be described by T-matrix the- 
ory, which treats the scattering of quasiparticles in a 
homogeneous (impurity averaged) medium off a single 
impurity in a self-consistent fashion. The resulting con- 
stant iV(0) then arises from a constant broadening 7 (de- 
fined above) of states near the d-wave nodes. In contrast, 
the low energy DOS in the BdG theory comes from new 
states arising out of hybridization of impurity resonances. 

We already see from Fig. 2 that at and beyond the 
critical concentration of the T-matrix approach, n.f„jp — 
0.08 for our choice of parameters, the BdG DOS does 
not approach the non-disordered value A^o- This raises 
the questions: does SC persist beyond nf^^, and if so, 
how? To address these issues we calculate the super- 
fiuid stiffness using the linear response result: -D^/tt = 
{^kx) — ^xxilx = 0, (7j, ^ 0, w = 0). The diamagnetic 
term (— fcx), is one-half (in 2D) the kinetic energy (— /C), 
and the paramagnetic term K^x, is the long wavelength 
limit of the transverse current-current correlation aver- 
aged over disorder realizations. 
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FIG. 4. r = (a) off-diagonal long order parameter 
and (b) superfiuid stiffness, as a function of concentration 
of unitary scatterers, obtained by the BdG method. Note 
that d-wave superconductivity is much more robust than the 
T-matrix prediction. Parameters used are J — U = 1.15t 
and (n) = 0.875, with A^o = 0.21 and Ds,o = 0.80, on an 
A'^ = 24 X 24 system, averaged over 15 disorder realizations. 

We see from Fig. 4b that the superfiuid stiffness Dg 
is much larger than the T-matrix result, consistent with 



Ref. PI, and does not vanish up to ni,„p = 0.12 which 
is 50% larger than nf^p within the T-matrix approxima- 
tion. (We did not go to higher impurity concentrations 
because of the increase in computational time to reach 
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self-consistency.) In any case, we expect that once Ds is 
sufficiently small, phase fluctuations neglected within the 
BdG mean field approach will drive the transition to the 
non-superconducting state this is left for a future 
investigation. Here we wish to gain insight into how the 
system manages to exhibit Ds > 0, even when T-matrix 
theory predicts it to be non-superconducting. 

One way to think about this is to correlate Dg and 
N{uj). A smaller DOS for low- lying excitations in the 
BdG approach implies fewer "normal fluid" excitations 
and hence a larger superfluid density compared to the 
T-matrix approximation. A complementary approach, 
which we find very illuminating, relates the Ds to the 
inhomogeneous pairing amplitude A{i) in the disordered 
ground state, shown in the left panels of Fig. 3. No- 
tice that the d-wave pairing amplitude is suppressed in 
the vicinity of an impurity on the scale of the coherence 
length ^0 which is 3 to 4 lattice units. (In addition, a 
small extended s-wave component, not shown, also de- 
velops nearby) . The regions of suppressed pairing ampli- 
tude give the appearance of "swiss cheese" Q at finite 
"imp in Fig. 3. 

In the T-matrix approach the order parameter is forced 
to be spatially uniform and it vanishes for riimp > '^fmp- 
However, by allowing the pairing amplitude to vary on 
the scale of ^o, in response to the impurity potential, the 
BdG solution permits a non-vanishing order parameter 
A which is larger than that obtained within T-matrix 
theory for all niaip, see Fig. 4a. (A is formally defined 
in terms of the long distance behavior of the appropriate 
reduced two-particle density matrix) . We note that both 
A/Aq and Ds/Dsfi are linear functions of n-i^p^Q for a 
substantial range of impurity concentration. 

To qualitatively understand the superfluid stiffness Ds 
consider applying an external phase twist to the inho- 
mogeneous ground state. Despite the fact that at large 
iT-imp there are large regions where the amplitude van- 
ishes, there are still paths that permit phase information 
to be conveyed from one edge of the system to the other, 
thus leading to a non- vanishing Ds. Thus the spatial 
inhomogeneity of the pairing amplitude, which is partic- 
ularly important in short coherence length superconduc- 
tors, is crucial in understanding the relative insensitivity 
of the system to unitary impurities, in that the order 
parameter and superfluid stiffness are much larger than 
one might have guessed from the T-matrix approxima- 
tion. This lack of sensitivity of the high Tc cuprates to 
disorder has been seen in numerous experiments 18|. 

Despite the quantitative results on finite systems and 
their detailed qualitative understanding, many questions 
remain open. The first one relates to Tc suppression. 
While it is easy to calculate the "mean field Tc" , a more 
reliable estimate should include the effect of both phase 
fluctuations and quasiparticles. Another important ques- 
tion is thermal transport [|l9| in the SC state. Why does 
it not reflect the low energy structure of the DOS and why 
is it consistent with_the universal behavior predicted by 
when the superfluid density jl^ ] 



shows deviations from it. A full understanding of the 
asymptotic DOS of the low-energy excitations, their lo- 
calization properties and the study of SC state transport 
on a network of hybridized resonances are all topics for 
future research. 
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